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Abstract 

The nonequilibrium magnetic domain structure of growing ultrathin ferromagnetic 
films with a realistic atomic structure is studied as a function of coverage and tem- 
perature. We apply a kinetic Monte Carlo method to a micromagnetic model de- 
scribing the transition from superparamagnetic islands at low coverages to a closed 
ferromagnetic film. The magnetic relaxation and the island growth happen simulta- 
neously. Near the percolation threshold a metastable magnetic domain structure is 
obtained with an average domain area ranging between the area of individual mag- 
netic islands and the area of the large domains observed for thicker ferromagnetic 
films. We conclude that this micro-domain structure is controlled and stabilized by 
the nonuniform atomic nanostructure of the ultrathin film, causing a random inter- 
action between magnetic islands with varying sizes and shapes. The average domain 
area and domain roughness are determined. A maximum of the domain area and a 
minimum of the domain roughness are obtained as a function of the temperature. 

Key words: magnetic domains, nanostructured film, nonequilibrium states, 
random magnet, ultrathin film growth, kinetic Monte Carlo method 
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1 Introduction 



The investigation of the magnetic domain formation, in particular in nano- 
structured ultrathin films, is an active field of current research. The influence 
of the atomic morphology on the magnetic properties is known to be especially 
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strong during the initial states of the thin film growth. A small variation of 
the preparation conditions may change the corresponding magnetic structure 
markedly. This has been shown by the recent progress of highly resolving imag- 
ing techniques, which allows for the investigation of the atomic structure as 
well as of the magnetic domain structure of these nanostructured systems in 
greater detail [1-6]. For example, a strong dependence of the domain struc- 
ture on the film morphology has been observed for the Co/Au(lll) thin film 
system yielding a smaller average domain area for rougher Co films [1,2]. In 
addition, we emphasize that the consideration of nonequilibrium states for the 
investigation of the magnetic domain structure is of particular importance. 

In this contribution we study the interplay of the atomic and magnetic struc- 
ture theoretically. In particular we simulate the spatial and temporal magnetic 
domain formation during the simultaneous growth of ultrathin ferromagnetic 
films. To our knowledge this problem has not been investigated previously. 
The domain structure is calculated within a kinetic Monte Carlo (KMC) sim- 
ulation applied to a micromagnetic model considering exchange coupling and 
uniaxial lattice anisotropy. The average domain area S and the average domain 
roughness R of the growing thin film are calculated as functions of coverage 
© and temperature T. The domain roughness is defined as the edge-to-area 
ratio of a domain. 

In general, for low coverages the growing thin film mostly consists of isolated 
superparamagnetic islands with a domain area of the order of the island area. 
On the other hand, at higher coverages the coagulation of islands results in 
a closed ferromagnetic film which allows for the formation of large domains 
with sizes of the order of about 1 yum, as usually observed [1]. The crossover 
between these two extremal domain size ranges is still not well investigated. 
Also, the quantities determining this transition are not known. 

For the simulation of the magnetic properties of the growing thin film we 
consider the following important features: 

i) We take into account an irregular, nanostructured film morphology which is 
usually present in a realistic thin film system. This is performed with the help 
of a simple growth model, allowing for the consideration of a size range between 
the atomic and the |Um-scale. The main quantities determining the growth 
(island density and -arrangement, growth mode, etc.) are chosen in accordance 
with experiments. Due to the strong direct exchange coupling each island is 
assumed to be single-domain, its magnetization rotates thus coherently [7] . 

ii) In the coverage range near the percolation threshold 0p of the thin film 
the interplay between the (irregular) atomic morphology and the magnetic do- 
main structure is expected to be strongest. Near 0p the domain formation is 
governed by 'weak links' between the magnetic islands, while the island-island 
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interaction tries to reach a ferromagnetic (single-domain) state of the film. The 
islands of a homogeneous, periodic island ensemble merge at the same time, 
thus the crossover from small to large domains takes place simultaneously all 
over the growing film. The characteristic times to cross the energy barriers, de- 
pending on the anisotropy and the island size, are the same for all islands. For 
an irregular nanostructure, however, the interactions between different island 
pairs will be strongly nonuniform (random magnet), causing various energy 
barriers for the magnetic rotation of different islands, and various inherent 
characteristic times [8,9]. A pre-formation of domains with intermediate areas 
covering several neighbouring islands is expected. This property infers much 
longer magnetic relaxation times to reach the thermal equilibrium than for a 
periodic island array. In other words, such a random ferromagnetic system can 
be more easily trapped into metastable states with a frozen-in multi-domain 
structure, and with an average domain size in between the above mentioned 
extremal sizes. Here the consideration of a realistic irregular atomic structure 
is especially important. 

iii) A finite temperature is taken into account in order to investigate the do- 
main formation of a realistic thin film system. The temperature affects strongly 
the probabilities of the island magnetizations to overcome the energy barriers 
resulting in a long range magnetic order depending on time and temperature. 
In addition, also the internal (short range) magnetic order of a single island 
due to the finite exchange coupling is considered, leading to temperature de- 
pendend island-island interactions and lattice anisotropics. 

iv) Furthermore, we emphasize that for the investigation of such an ensemble 
of magnetic islands the consideration of states far from thermal equilibrium is 
very important. The long relaxation times common to a random magnet may 
be much larger than typical detecting times. We expect in particular that 
in the coverage range around ©p the domain structure with an intermediate 
average domain area may be observable within the measuring time or the 
simulation time. 

v) Most importantly, the equilibrium state also changes permanently, if during 
the magnetic relaxation the thin film grows steadily. Thus, within our simu- 
lation of the magnetic domain stucture we take simultaneously into account 
the temporal variation of the atomic structure during the thin film growth. 
This results consequently in a temporal variation of the magnetic interactions 
determined by the growth velocity. 

The formation of large magnetic domains is expected to be facilitated by an 

increasing temperature, since then energy barriers can be more easily sur- 
mounted. This will result in a larger average domain area S after a given 
observation time, or a faster formation of larger domains. On the other hand, 
fluctuations at elevated temperatures will destroy large domains yielding a 
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smaller average domain area. Thus, a maximum of S{T) as a function of tem- 
perature during domain formation is expected. With the same argument we 
expect a minimum for the average domain roughness R{T). 

Previous simulations using similar micromagnetic models have been performed 
mainly for the investigation of the domain structure during the magnetiza- 
tion reversal of thicker ferromagnetic films with a granular structure [10]. In 
these studies the systems have been divided into periodic arrays of identical 
magnetic cells with a nonuniform distribution of the magnetic interactions be- 
tween the cells. In contrast, within our model we consider an irregular atomic 
structure. Different anisotropics and magnetic island-island interactions are 
caused by the island size dispersion and different common surface areas be- 
tween neighboring islands. To our knowledge such a realistic inhomogeneous 
atomic structure of an ultrathin film in the early states of film growth for 
the investigation of the spatial domain formation has not been considered 
previously. 

In addition, theoretical investigations often consider equilibrium magnetic 
structures [11]. As mentioned, the relaxation times to reach the equilibrium 
states may be very long, so that in experiments metastable multi-domain struc- 
tures are usually observed. In the present study we will simulate such domain 
patterns by calculating the temporal variation of the magnetic structure far 
from equilibrium. Since analytical approaches are not available for the study 
of the domain formation in a large system far from thermal equilibrium, we 
are forced to apply a kinetic (time resolved) Monte Carlo simulation. 

This work is structured as follows. In Sec. 2 we outline our model. A sim- 
ple growth model is applied which allows for the consideration of large sys- 
tems. Magnetic interactions are taken into account by a micromagnetic model, 
the magnetic relaxation is determined within a kinetic Monte Carlo method. 
Results for the average domain area and domain roughness as functions of 
coverage and temperature are presented in Sec. 3. A conclusion is given in 
Sec. 4. 



2 MODEL 

Eden Growth Model: We apply a Monte Carlo (MC) method [12] to sim- 
ulate the molecular beam epitaxial (MBE) growth of a magnetic ultrathin 
film on a nonmagnetic substrate using a simple solid-on-solid growth model. 
Within this so-called Eden model [12,13] each atom is randomly attached to 
already existing islands, and stays immobile afterwards. To take into account 
different growth modes we assume furthermore that the adatoms are placed 
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on lattice sites i with varying probabilities 

p{qi,z) (xexp{-A{z),/qi) . (1) 

Here z refers to the layer index. The square-root dependence on the local co- 
ordination numbers qi has been obtained to be approximately valid for metal 
surfaces [14]. By using layer dependent binding parameters A{z) we are able 
to simulate with simple means different growth modes such as an island-type 
growth mode (three-dimensional islands) or a layer-by-layer growth mode. 
Also various surface faces and island densities and -arrangements can be 
considered in accordance with experiments, for example islands arranged in 
chains. This modified Eden growth model is valid for the case of a fast surface 
diffusion and a moderate step diffusion for which the mobility of the atoms is 
large enough to probe the different atomic positions of the island edges. Since 
the very time consuming calculation of the atomic diffusion is avoided, we are 
able to consider a large number of nonequivalent sites (~ 10®). The influence 
of the temperature on the thin film morphology is not considered within our 
study [15]. 

In the present calculations we assume a 500 x 500 atomic and magnetic unit 
cell (2.5 • 10^ sites) on a fcc(OOl) lattice with periodic boundary conditions. 
We choose as an example the bilayer growth mode of the first two atomic 
layers as observed for the Co/Cu(001) thin film system [16], using always 
the ratio of the binding parameters A{l)/A{2) = 0.989. In accordance with 
experiment the island density is put equal to 0.005 per lattice site, refering 
to 1250 elementary islands in the unit cell. The simulation is started with a 
random distribution of occupied sites with minimal mutual distance rmin = 10 
lattice constants which serve as seeds. Each additional atom is placed on a 
perimeter site with probability p{qi, z). We stop the growth after deposition of 
two magnetic layers. The full coverage represents a system with two completely 
filled layers. 

Magnetic Structure: For a given atomic structure during the growth of the 
thin film the corresponding magnetic properties are determined by performing 
a kinetic Monte Carlo simulation [12]. The following micromagnetic model for 
the total energy of a system of interacting magnetic islands is assumed: 

^ = -7 E ^iAT, ©) hjiQ) Sj - y: k,{t, e) N,ie) {stf , (2) 

where Si is a normalized vector characterizing the direction of the magnetiza- 
tion of the i-th island, T is the temperature, and the coverage of the film. 
Each magnetic island with iVj atoms is assumed to be single-domain (Stoner- 
Wohlfahrt particle [7]) with a single giant magnetic moment /li — Ni /Xat, A*at 
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the atomic magnetic moment. The island moments are subject to a uniaxial 
lattice anisotropy Ki per atomic spin. Due to this anisotropy we simplify our 
calculations by considering only two different directions along the easy axis for 
each island moment {Si — ±1). This corresponds to the moments directed 'up' 
and 'down' for a perpendicular film magnetization, or to the two orientations 
of the in-plane magnetization in case of a rectangular (110) layer, for exam- 
ple. For an isolated island these two states are separated by an energy barrier 
given by its total lattice anisotropy AEi = Ki Nj, which has to be surmounted 
during magnetic reversal [17]. During the thin film growth two islands i and 
j may merge. If the directions of their magnetizations are different, the two 
islands minimize their mutual exchange coupling by the creation of a magnetic 
domain wall with the energy E = '-fijLij. Here '-fij is the domain wall energy 
per bond and L^jiQ) is the common surface area of the islands in units of 
the lattice constant a-o [18]. For simplicity the magnetic dipole coupling is not 
considered here. 

We emphasize that the energy expression Eq. (2), representing a system of 
individual magnetic islands with varying interactions jijLij, is a good ap- 
proximation for nanostructured systems. By application of this model we can 
describe isolated superparamagnetic islands at low coverages (Ly(©) = 0), 
as well as a strongly connected ferromagnetic film at high coverages, and the 
transition between these two extremes during the thin film growth. However, 
in the case of a smooth film (in the present study for ^ 1.8) the mainte- 
nance of irregularly shaped magnetic islands with nonuniform couplings is an 
unphysical discretization of the system. Thus in such a coverage range our 
model is not valid. 



Furthermore, the decreasing internal magnetic order of the islands with in- 
creasing temperature results in temperature dependend island interactions 
^ij{T) and effective anisotropy coefficients Ki{T). We treat the atomic mag- 
netic moments as S = 1 quantum spins. Within a mean field theory the relative 
internal island magnetization mj(T, 0) is given by the Brillouin function 



mi(r,e) = >SSi(a;,) = ^coth(^^a;i) -^cothQx,) , (3) 

_ ^.(e) Jm,(r,e) 

Here, ^^(0) is the average coordination number of island i which depends on 
the coverage © of the growing thin film, /cb is the Boltzmann constant. The 

ability of the anisotropy to maintain a certain direction of the magnetization 
decreases due to thermal agitation. Thus, a decreasing mj(T, 0) causes also 
a decreasing Ki(T). Within a first order thermodynamic pertubation theory 
[19] the effective anisotropy Ki{T) for 5" = 1 is given by 
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K 



4 — - coth [ -Xi] coth ( - ) + - coth^ ( - Xi 



(5) 



Note that mi{T) 1 and fi{T) 1 for T — 0. Following the general 
derivation of the domain wall width [20], the influence of the finite temperature 
on the interaction between the islands is considered as 



(6) 



where /ij(T, 9) is the average of /i(T, 9) and fj{T, 9), see Eq. (5), of islands 
i and j. 

The temporal development of the magnetic arrangement of the island ensemble 
is determined as follows. Thermal activation and (surface) interaction energy 
between neighboring islands may cause reversals of the island magnetizations 
between the two states = ± 1. This relaxational behavior during the film 
growth is calculated with the help of the KMC method [12]. From Eq. (2) the 
energy of island i as function of the angle (p between Si and the easy axis is 
given by 

ei(0) = Ei/Ki Ni = -2hi cos - cos^ , (7) 



with hi — SiJ2jlij Lij Sj/AKiNi. Here we have made use of the condition 
Sj — ±1. By analysing Eq. (7) two cases have to be distinguished. 

First for < 1 the two states Si = ±1 {(p = 0,7t) both represent energy min- 
ima. These are separated by an energy barrier governed by the competing vol- 
ume energy (lattice anisotropy) and surface energy contributions. The barriers 
for the transitions {Si = +1) ^ {Si = — 1) are given by AEi = Ni Ki{hi ± 1)^. 
The fiip rate Fj of island i to overcome A£^j is calculated by use of the 
Arrhenius-type ansatz [21] 



where the prefactor Fq is treated as a constant 'attempt frequency', determin- 
ing the time unit of the magnetic relaxation. 

Secondly, for \hi\ > 1 one of the states Si — ±1 refers to an energy maximum 
and the other to a minimum. Then AEi = ±4A^j Ki hi is the energy difference 
between the two directions of the island magnetization and is governed only 
by the surface energy. This case is treated with the usual Metropolis algorithm 
[12,22], using the same prefactor Fo as in Eq. (8). 
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The growing thin film is characterized by a large amoTint of noncquivalent lat- 
tice sites, corresponding to a large number of different interaction parameters. 
Since little is known about these values, we use in our simulation averaged 
quantities for 7, J, and K which are fixed as follows, using as an example 
the Co/Cu(001) thin film system. The domain wall energy 7 is adjusted to 
give the observed Curie temperature of the ferromagnetic long range order of 
Tc = 350 K of a Co/Cu(001) film with two monolayers (ML) [23]. We obtain 
7 = 5.8 meV per bond [18]. In a nanostructured film consisting of connected 
islands the internal ordering temperature of the islands T^* is larger than Tc- 
Due to the lack of measurements for T™*, the exchange couphng J is chosen 
to yield T™* = 450 K at a coverage of 2 ML, refering to a coordination number 
z = 8. This results in J = 7.3 meV per nearest neighbor bond. The uniaxial 
anisotropy constant is set equal to K = 0.2 meV per atom, which is a typical 
value for ultrathin films [24] . For comparison, with these values the width of 

an undisturbed Bloch wall in a bulk ferromagnet is given by wb = tt J/2K fa 

13 Bio, and the corresponding wall energy 73 = 4^ JK/2 = 3.4 meV per bond 
[20]. The 'attempt frequency' L^ is set equal to = 5.0 • 10^ sec"^ [21]. The 
growth velocity of the film is adjusted to 1 ML/ 100 sec. 

The simulation procedure is performed as follows: first the island seeds with 
random spin directions are distributed over the unit cell with periodic bound- 
ary conditions. Then for a given temperature the thin film grows up to a 
coverage of 2 ML, which is divided into 400 growth steps. After each growth 
step the quantities Ni{Q), Lij{Q), '~fij{T,Q), and Ki{T,Q) are calculated. By 
applying 100 Monte Carlo steps (MCS) per island moment of randomly cho- 
sen islands the magnetic structure of the island ensemble is allowed to relax 
towards its thermal equilibrium. Then the growth procedure is repeated. 

The magnetic domains are determined by identifying connected islands with 
a parallel magnetization (Hoshen-Kopelman algorithm [25]). From this pro- 
cedure the number of domains as well as their average area S and average 
roughness R are calculated. Only the first magnetic layer is used for this anal- 
ysis. To indicate the influence of the atomic structure on the magnetic domain 
properties, we calculate also the average area of the elementary islands, as well 
as the average area of connected islands, consisting of coagulated elementary 
islands. 

The domain roughness Ri is defined as the ratio Ri = dSi/Si of the number 
of atomic magnetic moments dSi at the outer edge of the domain and its 
total number S^. Since for a domain with an increasing area its edge-to-area 
ratio is always decreasing, we consider here rather the average relative domain 
roughness r = R/Ro with respect to the minimal roughness (smallest possible 
edge) RI oc S^^^"^ of a circularly shaped domain with the same area S^. 
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3 Results and Discussion 



Using the Eden growth model, the KMC procedure, and the interaction param- 
eters as described in the preceding section we present results for the magnetic 
domain structure of a growing ultrathin film. In particular, we have deter- 
mined the average domain area and domain roughness as functions of the 
coverage and the temperature. The results are always averaged over 20 runs 
for systems with same global parameters determining the thin film growth and 
the magnetic properties. 

Examples of the resulting atomic structure as obtained from our growth pro- 
cedure are shown in Fig. 1, using the above described bilayer growth mode. 
For a better visuahzation of the film morphology a unit cell with 300 x 300 lat- 
tice sites is chosen here. Snapshots of three different coverages are depicted, 
e = 0.5 ML, 1.0 ML, and 1.5 ML. The coverage = 2.0 ML corresponds 
to a smooth magnetic film with two closed layers (not shown). The result- 
ing atomic structures are similar to those observed experimentally for the 
Co/Cu(001) system [16]. 

In Fig. 2 we show snapshots of resulting magnetic domain structures for dif- 
ferent coverages and for the temperature T = 100 K. As can be seen, for 
small coverages the domain resembles the island For large cover- 

ages 0^2 ML the thin film is almost closed and consequently the magnetic 
domains are very large. Very interestingly, for coverages near and above the 
percolation threshold 0p --^ 0.9 ML the domains assume an intermediate size, 
covering several neighboring islands. We claim that this micro- domain struc- 
ture is controlled and stabilized by the nonuniform atomic nanostructure of 
the growing ultrathin film. 

From the atomic and magnetic structures as depicted in Figs. 1 and 2 we 
analyse now the average domain area and domain roughness as functions of 
coverage and temperature. The average domain area S in units of lattice sites 
is given in Fig. 3(a,b) as functions of the coverage and the temperature T. 
The average elementary island area -S'isiand and the average area of connected 
islands, which serve as lower and upper limits of S, are also shown in Fig. 3(a). 
For coverages well below the percolation threshold 0p the average domain area 
is of the order of S'isiand) see Fig. 2. In this case the ultrathin film consists of 
magnetically almost isolated islands, the exchange coupling between neigh- 
boring islands has no large infiuence. If no additional long range magnetic 
interactions such as the dipole coupling or the indirect exchange (RKKY-) 
coupling are important, the system of isolated magnetic islands refers to a 
supcrparamagnet. Due to their small size and their superparamagnetic be- 
haviour the magnetic domains in this coverage range might be hardly visible 
with experimental means. 
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With increasing coverage the islands start to merge and form large connected 
islands. Caused by the nonuniform island ensemble the island coalescence does 
not occur simultaneously in the system. In this coverage region the average do- 
main area increases considerably, but is still markedly affected by the atomic 
nanostructure of the thin film. We find an average domain area definitely 
larger than the average elementary island area, but still much smaller than 
the large domain areas for an almost closed ferromagnetic film. Due to the 
irregular atomic structure of the island ensemble, the area of the interfaces 
between neighboring islands differ considerably, resulting in strongly different 
surface contributions of the energy barriers. Such a distribution of the effective 
interactions within the island ensemble can be attributed for by a 'random ex- 
change energy' between neighboring island magnetizations. For such a random 
magnet metastable spin structures exhibit considerable temporal stabilities, 
characterized by large magnetic relaxation times [8,9]. Only after a long time 
the system may overcome the energy barriers and reach its equilibrium mag- 
netic state. We expect that this mechanism also explains the appearance and 
stability of the magnetic micro-domain structure above the percolation thresh- 
old of the thin film as obtained from our calculations, since here the irregular 
atomic morphology causes long magnetic relaxation times. 

In case of a uniform magnetic island structure, including also the case of 
a smooth ferromagnetic film, the corresponding relaxation times are much 
shorter and a micro-domain structure might not be observable. We have tested 
this assumption by simulating the domain structure also for a periodic array 
of identical magnetic islands with equal interactions between them, in con- 
trast to an irregular island system. As expected, close to 0p the domains 
increase fast to a large average size. Thus, we conclude that the frozcn-in 
micro-domain structure with an intermediate average domain area is caused 
by the nonuniform island ensemble. The experimentally observed small do- 
mains in ultrathin films with a rough morphology [1,2] are expected to be 
such a metastable, nonequilibrium magnetic structure. 

With further increase of the coverage a smooth magnetic thin film is obtained, 
and the domain pattern is mainly determined by competing magnetic inter- 
actions (exchange, anisotropies) . The domain area distribution is quite wide, 
few large domains coexist with a number of small domains. Due to the dis- 
appearance of these small domains, see Fig. 2(c,d), the average domain area 
becomes very large. Since the extension of these large domains reaches the 
size of our unit cell, we emphasize that in the coverage range of about 1.8 - 
2 ML our simulation of the domain structure is influenced by finite size effects. 
As discussed in the previous section, our model should not be applied for a 
smooth film, which is present in this coverage range. 

For the temperatures as considered in Fig. 3(a), for most coverages 9 a larger 
average domain area S{Q) is obtained for a smaller T. Investigating S{T) as 
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a function of the temperature during the domain formation, see Fig. 3(b), we 
obtain a maximum of S{T) in particular at low T. The reason is that for low 
temperatures the domain formation is hindered by energy barriers. Since with 
an increasing T the probability to overcome energy barriers becomes larger, 
an increasing average domain area is obtained. On the other hand, due to 
thermal agitation a large domain may disintegrate into domains with smaller 
sizes, resulting in a decreasing average domain area with increasing tempera- 
tures. While approaching the Curie temperature T — > Tc(2ML) = 350 K the 
average domain area drops considerably. With an increasing coverage the 
maximum of S{T) is shifted to larger temperatures, since the increasing aver- 
age magnetic energy between the islands is better suited to withstand thermal 
agitations. 

In Fig. 4(a,b) we present results for the average domain roughness as functions 
of the coverage © and the temperature T. The domain roughness is charac- 
terized by the edge-to-arca ratio of the domains. As described in Sec. 2, we 
present here the relative domain roughness r = R/ Ro with respect to the 
minimal roughness R^ of circularly shaped domains. Values f < 1 result from 
single-domain states of the simulated films. Similar as for the average domain 
area S{T), an increasing temperature facilitates the surmount of energy barri- 
ers, and the domains may assume a more compact average shape, i.e. arc less 
rough. In contrast, an even higher temperature will destroy ordered structures 
in particular at their surface. Consequently, a minimum of r as function of T is 
observed, see Fig. 4(b). Evidently, the relative roughness depends also on the 
average domain area and the connectivity of the system, as can be seen from 
Fig. 4(a). An increasing temperature tends to smoothen the domains with a 
large area at larger coverages 0. For lower coverages, however, the opposite 
behavior is observed: a larger roughness for an increasing temperature. This 
might be explained by the fact that in this coverage range the average do- 
main area and the connectivity between the islands are comparably small so 
that the disordering effect resulting from thermal fluctuations dominates the 
smoothening due to a facilitated surmount of energy barriers. 



4 Conclusion 

In this study we have investigated the magnetic domain structure of an ul- 
trathin film in the early states of film growth. Here the infiuence of the film 
morphology on the domain formation is especially strong. A model has been 
applied which is suited in particular for nanostructurcd systems, as well as 
for the consideration of nonequilibrium states. To our knowledge, the occur- 
rence and relative stability of the domain structure has not been considered 
previously in connection with a growing magnetic thin film [26] . 
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The average domain area and the average relative domain roughness have been 
calculated as functions of the coverage and the temperature. In particular, we 
find an interesting metastable micro-domain structure near the percolation 
threshold of the thin film. The average area of these domains ranges between 
the area of the elementary islands and the large domains of a smooth ferro- 
magnetic film. We conclude that this nonequilibrium micro-domain structure 
is controlled and stabilized by the nonuniform spatial (nano-) structure of 
the thin film. Such a mechanism could lead to the unusually small domains 
observed in rough ferromagnetic films [1,2]. Similar as in a random magnet 
a distribution of the magnetic couplings, as present in an irregular system, 
causes large magnetic relaxation times. For a thin film system with a peri- 
odic array of magnetic islands (positions, sizes, shapes, interactions), a micro- 
domain structure might not be observable, since large domains will evolve fast 
due to a short relaxation time. The average domain area is found to increase 
with increasing coverage, and to exhibit a maximum as function of the tem- 
perature. Similarly, according to our calculations the average roughness shows 
a maximum as function of the coverage and a minimum as function of the 
temperature. 

We have considered single-domain magnetic islands with a coUinear magne- 
tization during reversal (Stoner-Wohlfarth particles [7]). This approximation 

is valid if the island size is well below the domain wall width w oc \J J / K 
of a bulk ferromagnet [20], yielding an upper limit for the island size within 
our calculations. For larger sizes a noncoUinear island magnetization has to 
be taken into account [27]. Furthermore, as discussed in Sec. 2 the present 
model with nonuniform island interactions should not be applied for a smooth 
ferromagnetic film. Nevertheless, we expect that our approach is also suited 
for the investigation of thicker films with a considerable surface roughness. 

Also other thin film magnetic quantities will be affected strongly by an irregu- 
lar nanostructure and the corresponding large magnetic relaxation times. For 
instance, the remanent magnetization as measured e.g. by MOKE, and the 
magnetization reversal by an applied magnetic field may depend considerably 
on the nonuniform atomic structure [28]. Since in near future magnetic nanos- 
tructured systems with a defined lateral structure will be prepared in a more 
and more controlled manner [29], the corresponding magnetic properties can 
be calculated for these atomic structures. Vice versa, information about the 
atomic morphology can be obtained by analysing the magnetic properties. 

Some remarks concerning the improvement of our present model study are 
in order. First, we will not only consider the two directions (Sj ± 1 of the 
magnetic islands, but will allow for a continuous rotation of the magnetiza- 
tion. This will result in particular in a noncollinear magnetic arrangement of 
the island ensemble. The formation of larger magnetic domains may be facili- 
tated, resulting in a faster magnetic relaxation. Secondly, the magnetic dipole 



12 



coupling will be included in our model, leading possibly to an even more com- 
plex magnetic structure. Thirdly, we will calculate explicitely the relaxational 
behavior of the magnetization for different coverages and temperatures of a 
nanostructured thin film. Prom an analysis of the temporal behavior the re- 
laxation times can be extracted. Generally, also the growth velocity will affect 
the domain structure: for small velocities the islands have more time to align 
themselves and to create larger domains or a single-domain state. On the 
other hand, for a fast growth the islands become too fast too big to overcome 
the energy barriers, and the resulting domains are expected to be smaller 
[30]. An applied magnetic field will support the formation of a single-domain 
state. Also, we will calculate the susceptibility in response to a magnetic field, 
which will yield additional informations for the interplay between the atomic 
and the magnetic structure of growing films. Furthermore, the temperature 
dependence of the atomic growth can be included by applying appropriate 
growth mode parameters. In the present study we have used the same growth 
mode for all temperatures in order to concentrate solely on the temperature 
dependence of the magnetic properties. Finally, we will improve our KMC 
method by allowing for simultaneous fiips of connected islands (cluster-spin- 
fiip algorithm [12,31]). The consideration of such correlated jumps will yield 
a more realistic relaxational behavior. 
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Fig. 1. Simulation of the atomic structure of a thin film with a bilayer growth 
mode. The unit cell has 300 x 300 lattice constants and contains 450 islands (island 
density 0.005 islands per site). Snapshots of three different coverages are depicted, 
(a) e = 0.5 ML, (b) 9 = 1.0 ML, and (c) 6 = 1.5 ML. Black refers to the uncovered 
substrate, light gray to the first and dark gray to the second magnetic layer. 
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Fig. 2. Snapshots of magnetic domain structures at T = 100 K for different coverages 
(a) e = 0.5 ML, (b) 9 = 1.0 ML, (c) 6 = 1.5 ML, and (d) 6 = 2.0 ML. The unit 
cell has 500 x 500 lattice constants and contains 1250 islands (island density 0.005 
islands per site). Only the domain pattern of the first magnetic layer is shown. The 
two gray scales refer to the two magnetic directions, the uncovered substrate is 
black. 
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Fig. 3. Semi-logarithmic plots of the average domain area S as functions of (a) the 
coverage B and (b) the temperature T during growth. In (a) the average elementary- 
island area and the average area of connected islands are the lower and upper limits 
of S. A nanostructured thin film system with the atomic and magnetic properties as 
described in the text has been assumed. The results are obtained from 20 different 
runs using the same growth mode and magnetic parameters. The area of the unit 
cell is 2.5 • 10^ sites, the percolation threshold is close to @p ~ 0.9 ML. 
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Fig. 4. Average relative domain roughness r = R/Ro as functions of (a) the coverage 
and (b) the temperature T during growth. The domain roughness R is given 
relative to the minimal roughness Ro of a circularly shaped domain with the same 
area as described in Sec. 2. The results are obtained from 20 different runs using 
the same growth mode and magnetic parameters. Values r < 1 result from single 
domain samples. 
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(a) 0=0.5 ML (b) 0=1.0 ML 




